Instability of Magnons in Two-dimensional Antiferromagnet at High Magnetic Fields 
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Spin dynamics of the square lattice Heisenberg antiferromagnet, Ba2MnGe2 07, is studied by a 
combination of bulk measurements, neutron diffraction, and inelastic neutron scattering techniques. 
Easy plane type antiferromagnetic order is identified at T < 4.0 K. The exchange interactions are 
estimated as J\ — 27.8(3)/xeV and J2 = 1.0(1) /xeV, and the saturation field He is 9.75 T. Magnetic 
excitation measurements with high experimental resolution setup by triple axis neutron spectrometer 
reveals the instability of one magnon excitation in the field range of 0.7 He < H < 0.85ifc- 



Many phenomena in condensed matter science can be 
explained by using the concept of quasiparticle. For ex- 
ample antiferromagnetic order is a result of Bose con- 
densation of magnons 1 -^ and superfluidity is those of 
phonons.- The quasiparticles in condensed matters, how- 
ever, can be unstable and decay if allowed by conserva- 
tion laws. Such phenomena was initially predicted in su- 
perfluid Helium and was identified by a termination of the 
excitation at twice the energy of a roton^ The magnon 
version of the spectral termination was observed in a few 
quantum magnets, where the one-magnon branch crosses 
the lower boundary of the two-magnon continuum. 5,6 Re- 
cently the instability of magnons with a similar mecha- 
nism is predicted in the 2D square lattice Heisenberg an- 
tiferromagnets (SLHAF) in high magnetic field^ - — At 
zero field a two-magnon continuum spreads in the higher 
energy region for all wave vector q and there is no de- 
cay channel for one magnon. With increasing field the 
one-magnon branch moves higher around q ~ (ir n) and 
eventually overlaps with the continuum at a threshold 
field H* ~ 0.76i^c- The hybridization of one-magnon 
with two-magnon continuum induces instability of the 
one-magnon state. While most theoretical calculations 
including spin- wave ^ Monte Carlo, 8 and exact diag- 
onalization 10 are for the spin S = 1/2 case, the field 
induced magnon instability is universal to 2D and 3D 
AFMs with general spins. It is the change of the cur- 
vature of AF magnon dispersion from downward at low 
field to upward at higher field that qualitatively explains 
the kinematical instability of the magnon^ 

Many model compounds for SLHAF have been re- 
ported so farii The required specification for experimen- 
tal study of the spin dynamics in high field is moderate 
exchange energy and reasonably low saturation field that 
can be reached by superconducting magnet. From this 
one may naively infer that a compound with longer inter- 
magnetic- ion distance, e.g., Ba2CuGe2 07 with 6.0 Ai 2 - is 
favorable for experimental investigation. Unfortunately, 
however, the measured Hq for Ba2CuGe207 was 27 T— 
which is slightly higher than achievable in most supercon- 



ducting magnets. Hence we considered the presumably 
isostructural Mn 2+ (S = 5/2) compound Ba2MnGe207 
as a better experimental candidate. In this letter we 
report the experimental observation of a general nature 
of antiferromagnetic magnon, namely, the field induced 
magnon instability in 5=5/2 SLHAF Ba 2 MnGe 2 7 by 
combination of bulk susceptibility, high field magnetiza- 
tion, and neutron scattering techniques. 

A single crystal with 7 x 7 x 80 mm 3 was grown by 
the floating zone method. The crystal structure was 
determined by using a commercial x-ray diffractometer 
with area sensitive CCD. This confirmed that the sam- 
ple has a tetragonal structure with space group P\2\m 
and lattice parameters with a = 8.505 A and c = 5.528 A 
which is indeed isostructural to Ba2CuGe2 07. As shown 
in Fig. [T](a) and (b) Mn04 tetrahedra linked by GeC^ 
tetrahedra forms a square lattice in the crystallographic 
c plane. The Mn 2+ -Mn 2+ distance is 6.014 A which in- 
dicates that the exchange constant may be reasonably 
small. The square lattice is separated by Ba 2+ ions in 
the crystallographic c direction and the system is pre- 
sumed to be 2D-SLHAF. 

Magnetic susceptibilities, x's, along [0 1] and [1 1 0] 
were measured by a commercial SQUID magnetometer 
as shown in Fig. H{c). The applied field is 100 Oe. At 
T > 4.0 K anisotropy is quite small and almost negligi- 
ble. The gyromagnetic ratios of Mn 2+ ions are obtained 
g = 2.005 and 2.002 for H || [1 1 0] and [0 1], respec- 
tively by room temperature ESR measurement. Sub- 
stantially small anisotropy is observed due to the or- 
bital angular momentum of Mn 2+ ion for which the high 
spin state is zero. This indicates that spin exchange 
interactions are purely isotropic in Ba2MnGe207. The 
broad maximum at T ~ 6 K indicates short range AF 
fluctuation. The susceptibility is reasonably fitted by 
5 = 5/2 2D-SLHAF 14 with J = 26 fieV (solid curves). 
Anisotropic behavior at T < 4 K indicates a magnetic or- 
dering, x m [1 1 0] does not go to zero at T = K. While 
measuring the susceptibility in the c plane, i.e., a* — a* 
plane, by rotating the sample around the c axis we ob- 
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FIG. 1. (Color online) (a) Crystal structure of Ba2MnGe207 
projected onto the crystallographic a plane. The solid rect- 
angle indicates one unit cell. (b) Crystal structure pro- 
jected onto the c plane, (c) Magnetic susceptibility along 
if || [1 1 0] and [0 1]. Solid curves are theoretical calcu- 
lation for £ = 5/2 SLHAFM.^ Both curves cannot be dis- 
tinguished within the printable resolution, (d) Magnetization 
curve at T = 0.5 K. (e) Temperature dependence of the mag- 
netic Bragg peak at q = (1 0.5). (f) Magnetic structure of 
Ba2MnGe2 07- Magnetic unit cell is doubled in the c direc- 
tion. Moments are confined in the c plane. 

served almost no anisotropy. This indicates an easy-plane 
type AF order. High field magnetization along H \\ [110] 
at T = 0.5 K was measured up to 14 T by using the 
pulsed magnet installed at IMR Tohoku university. The 
magnetization saturates at H > 9.75T = Hq, which is 
consistent with that estimated from S = 5/2 SLHAF 
with J = 26 /aeV. A peak structure ascribed to spin-flop 
transition was observed in the field derivative of mag- 
netization at H ~ 0.12 T. We attribute this to a small 
anisotropy in the c plane. 

Neutron diffraction measurements were performed at 
HQR spectrometer in JRR-3M, Japan Atomic Energy 
Agency (JAEA). Magnetic Bragg peaks are observed at 
(ft (n + l)/2) (n is odd) below T < 4.0 K = T N as 
shown in Fig. Die). 18 magnetic Bragg reflections were 
collected in the a plane. We found that easy-plane type 
AF order as shown in Fig. [T] (f) is consistent with the 
observed Bragg intensities. The direction of spins in the 
c plane is not determined by the unpolarized neutron 
experiment. The magnitude of the moment is estimated 
MMn = 4.66(6)/i B . 
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FIG. 2. (Color online) (a) Constant q scans at (ft 0) at T = 
1.7 K and H — 0T. A bar in the bottom shows experimental 
resolution, (b) Dispersion relations in zero field. Right axis 
is scaled by the spin wave band energy SJiS. 

Inelastic neutron scattering experiments in zero field at 
T = 1.7 K were performed at HER spectrometer in JRR- 
3M in JAEA. Collimation setup was guide-open-80'-80'. 
A Be filter was installed in front of sample and fixed fi- 
nal energy mode with Ef — 3 meV was used. Constant 
q scans at q = (ft 0) are shown in Fig. Efa). The bar 
in the bottom indicates the energy resolution ] 15 i 16 Well 
defined resolution limited peaks are observed at hu < 0.6 
meV. The excitation energy at q = (1 0) that corre- 
sponds to two dimensional antiferromagnetic wave vector 
q 2B = (tt 7r) is 0.1 meV. This suggests a weak interplane 
coupling along the c direction. After a series of scans 
in three directions the dispersion relation is obtained in 
Fig. [2Kb). Enhanced dispersion with the band energy of 
about 0.55 meV is observed in the c plane. Small dis- 
persion in the c* direction is identified. The obtained 
data is analyzed by classical spin wave theory where the 
exchange interactions J\ and J2 in Fig. [TJa) and (b) are 
considered in the Heisenberg Hamiltonian, 

H = J 1 Y,Si'S j + J 2 Y,Sk-S l . (1) 

(i,3) 0M> 

Here the sums are taken for pairs of spins. Obtained pa- 
rameters are J\ = 27.8(3)/ieV and J2 = 1.0(1) fieV. The 
value of J\ is consistent with that estimated by suscep- 
tibility measurement. Thus Ba2MnGe207 is confirmed 
as a quasi 2D SLHAFM with reasonable energy scale for 
observing the spin dynamics in high field. 

We performed inelastic neutron scattering with wide 
resolution for the initial study of high field spin dy- 
namics. The data were collected at LTAS spectrometer 
in JRR-3M in JAEA. Collimation setup was guide-80'- 
open-open. A Be filter was installed in front of sample, 
horizontal focusing analyzer and the Ef = 2.6 meV mode 
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FIG. 3. (Color online) False color plots of a series of constant 
q scans under magnetic fields with wide energy resolution. 
Solid curves are classical spin wave calculation. 



was used. The field was applied in the a direction and the 
temperature of the sample was T = 0.25 K. False color 
plots of constant q scans at q = (h 0) are shown in 
Fig. [3l With increasing field the energy of the main spin 
wave, mode 1, at q = (2 0) corresponding to q 2B = (0 0) 
increases and the band energy reaches about 1.1 meV at 
H = 10 T (~ Hq). The dispersion has a ferromagnetic 
profile with reduced unit cell in the a plane at H = Hq- 
In the intermediate field another spin wave, mode 2, is 
induced due to doubling of the magnetic unit cell in the c 
direction. At Hq the doubling is dissolved and the mode 
2 disappears. 

Solid curves are dispersions calculated by classical spin 
wave theory. The ground state of canted AF order with 
tilt angle from easy plane 6 = sin -1 (H / Hq) is assumed. 
Dispersion for mode 1 is obtained that 



■j(q))(j(0) -cos26 j(q)) (2) 



where j(q) is Fourier transformation of J' s. The parame- 
ters J\ and J2 estimated by zero field data are used. The 
calculation reasonably reproduces the data in all fields. 
In this wide experimental resolution no broadening was 
observed. 

To obtain the precise peak profile of the magnetic exci- 
tation we performed the measurements with high energy 
resolution. The data were collected at V2/FLEX spec- 
trometer in Helmholtz Zentrum Berlin. We put tight 
collimations, guide-60'-60'-60', with Ey fixed at 2.6 meV 
which gives the best resolution for practical use of 3-axis 
spectrometer. A magnetic field was applied in the c di- 
rection. Constant q scans at q = (h 0) with H = 8 
T and T = 0.06 K are shown in Fig. H At h < 1.3 
and h > 1.8 the excitations are within the calculated 
experimental resolutio n 15 ^ 6 shown by the shaded area. 
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FIG. 4. (Color online) Constant q scans at H 
areas indicates energy resolution. 
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Meanwhile at 1.4 < h < 1.7 we clearly observe a broad- 
ening of the peaks. The shoulder structures at h = 1.4 
and 1.6 are a contribution from another mode. The en- 
ergy widths are separately obtained by Gaussian fittings 
and plotted as function of h for H = 6, 8, and 10 T in 
Fig. [5fa). Solid curves are the experimental resolution 
at the peak energies. At H = 6 T and 10 T, the widths 
are resolution limited for all h. Significant broadening 
is identified at 1.4 < h < 1.7 at H = 8 T. In the field 
dependence of the width at q = (1.7 0) in Fig. [5jb), 
the broadening is observed at7T<i^<9T. 

The field induced instability of the magnon is generally 
predicted in 2D- and 3D-AFM. The necessary condition 
for the instability is that the one-magnon branch enters 
into the two-magnon continuum. Energy and momentum 
conservation laws for one-magnon and two-magnon arei 
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Here Q is antiferromagnetic vector. We numerically ob- 
tained the lower boundary of the 2-magnon continuum, 
^qX^' ky usm g Eq. j2j)-((lj) and the obtained parameters 
Ji 'and J 2 . At low field the condition c^mag > uj^z 
holds for any q and there is no decay channel. With 
increasing field Uq mag goes higher and the decay channel 
comes to open at H = H* . In Fig.[5tc) at H = 8 T we see 
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(red solid curve) is lower than cj^ mag (blue 



4 




(b) 0.16 
0.14 
?0.12 
| 0.10 
i 008 

Ll_ 

006 
004 



" C- q= (1.70 0) 
'_ — Resolution 

T=0.06 K Q 


5 


5 




- Q 

-o — 






D" 


i.i. 




i 


i 



0.2 ^ 



Field (I) 




FIG. 5. (Color online) (a) Energy width of the constant q 
scans at different fields. Solid curves indicate the experimen- 
tal resolution at the peak energies. In the right axis the en- 
ergy is scaled by spin wave band width at zero field SJ±S. 
(b) Field dependence of energy width of constant q scans at 
q = (1.7 0). (c) Blue curve is one-magnon energy and red 
curve is the lower boundary of calculated two magnon contin- 
uum. Symbols are the subtraction of the latter by the former. 
Dotted arrow indicates the decay channel allowed by conser- 
vation laws. 
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solid curve) in intermediate q. For clarity oo q m ^ 
is plotted by filled circles. It is identified that a decay 
channel opens at 1.22 < h < 1.70. This h range is totally 
consistent with h where the broadening of one-magnon is 
experimentally observed in Fig.[5fa). Hence the observed 
broadening of the peak profiles is ascribed to the field in- 
duced instability of magnon that has been theoretically 
prediced^"— 



For convenience the energy scaled by spin wave band 
width in zero field 8J\S is shown in the right axes in 
Fig. Efb), [5fa) and (b). The most enhanced broadening 
is about 25% of 8 Ji5 at H = 8 T - 0.8H C . The magnon 
spectrum of the theoretical calculation for 5=1/2 case is 
shown in Fig. 3 of Ref0 At H = O.S5H c a one-magnon 
centered sharp peak with broad two-magnon continuum 
is predicted. At H = 0.9Hq one-magnon is mostly de- 
cayed in wide q and the energy width is about 35% of 
8J\S. A similar result is also shown in e.g., Ref,8. Note 
that the energy in the theoretical papers is scaled only by 
J. The line width of Ba2MnGe207 is somehow compa- 
rable to 5 = 1/2 calculations. This is rather surprising 
because 5 = 5/2 spin is classical and hybridization of 
one- and two-magnons is expected to be small. Indeed 
in a short report^ 7 - the width for 5 = 5/2 SLHAF is 
calculated - J/10, only 0.5% of 8Ji5 at its best. The 
enhanced instability in Ba2MnGe2 07 might be related 
to the crossing of mode 1 and mode 2 that occurs in the 
same q regime in Fig. [3j 

At H > Hq with increasing field the decay chan- 
nel allowed by conservation law eqs.(j3]) and (j4j) be- 
comes larger in the q space. Consistently calculations 
for the 5 =1/2 case shows that the one magnon de- 
cays almost everywhere at H = 0.9 Hq£ It is only for 
H > 0.99i^c that the sharp one magnon excitation recov- 
ers due to the total suppression of quantum fluctuations. 
In Ba2MnGe2 07 the one magnon recovers much faster as 
shown in Fig. [Sfb). Already at H = 0.9Hc the peak is 
within our experimental resolution. The restricted field 
range is presumably because the 5 = 5/2 spin is more 
classical than 5 = 1/2. 

To conclude we observed field induced one-magnon in- 
stability in 5=5/2 2D-SLHAF Ba 2 MnGe 2 7 by inelastic 
neutron scattering. The instability field region is identi- 
fied as 0.7H C <H< 0.85# c at q = (1.7 0). The width 
of the broadened excitation is comparable to theoretically 
calculated 5 = 1/2 case. Advanced theoretical study for 
classical spin would be important. Magnon instability is 
now believed to be ubiquitous in various types of mag- 
net including triangular and kagome lattices. 18 Further 
experimental study in new compounds would be interest- 
ing. 
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